Quasiscarred modes and their branching behavior at an exceptional point 
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We study quasiscarring phenomenon and mode branching at an exceptional point (EP) in typically 
deformed microcavities. It is shown that quasiscarred (QS) modes are dominant in some mode group 
and their pattern can be understood by short-time ray dynamics near the critical line. As cavity 
deformation increases, high-Q and low-Q QS modes are branching in an opposite way, at an EP, into 
two robust mode types showing QS and diamond patterns, respectively. Similar branching behavior 
can be also found at another EP appearing at a higher deformation. This branching behavior of QS 
modes has its origin on the fact that an EP is a square-root branch point. 

PACS numbers: 05.45.Mt, 42.55.Sa 



Open quantum/ wave systems often exhibit interest- 
ing phenomena that have no correspondence in closed 
systems. Openness effects of the systems are important 
to understand the physics of them. Quasiscarring phe- 
nomenon is an example found in a spiral-shaped optical 
microcavity [l^ . Unlike the scarred mode [3, ^ , the qua- 
siscarred (QS) mode shows amplitude enhancement along 
a closed path, not a periodic orbit, due to the openness 
effects such as Goos-Hachen (G-H) shift [1,[5[ and Fresnel 
filtering (FF) effect Although QS modes have been 

found only in the spiral-shaped microcavity, the openness 
effects necessary for quasiscarring are so generic that QS 
modes should be found in other deformed microcavities. 

On the other hand, even for integrable systems, the 
openness explains the localized mode pattern along peri- 
odic orbits [8|. It is also noted that the openness effects 
responsible for quasiscarring still exist even in the regular 
regime. Therefore, it is natural to expect the existence 
of QS modes in a slightly deformed microcavity showing 
almost regular ray dynamics. 

One of inherent features of open quantum/wave sys- 
tems is the non-trivial energy surface near a degenerate 
point, called exceptional point (EP) [9|,[iq|7 at which two 
eigenvalues and eigenfunctions coalesce. The interest- 
ing physics related to the EP has been studied in various 
open systems lll-[l3t , and recently an EP has been found 
in optical microcavities [l3-[l6|. The EP is a square- 
root branch point mathematically, so neighboring two 
modes can become quite different, as belonging to differ- 
ent branches, on a passage by an EP [14]. 

In this paper, we study QS modes and their branch- 
ing behavior at an EP in typically deformed microcav- 
ities, quadrupole and stadium deformations. The QS 
pattern can be understood by short-time ray dynamics 
near the critical line, in both chaotic and almost regu- 
lar regimes. It is also shown that the whispering-gallery 
modes (WGMs) localized around the critical line are con- 
tinued into either the QS branch or the diamond branch 
as the deformation increases. A degenerate point (EP) 
of high-Q and low-Q QS modes is found to play a crucial 




FIG. 1: (Color online) (a) A QS mode in a chaotic micro- 
cavity. {A2 mode at s = 0.12) (b) Husimi function of the 
QS mode shown in (a), (c) A high-Q QS mode in an almost 
regular microcavity. {A2 mode at e = 0.06) (d) Husimi func- 
tion of the QS mode shown in (c). In (b) and (d), the PSOSs 
are shown in background and the diamonds and rectangles 
represent diamond and rectangle periodic orbits. 



role in the mode branching. 

Although the quasiscarring has been studied only in 
a spiral boundary shape [l|, [l3, [isl , we can expect that 
QS modes, formed by openness effects, would exist in 
other deformed microcavities. We thus take a typical 
deformation, quadrupolar deformation, given as r(0) = 
R{1 -\- £ cos 20) in polar coordinate with R = 1. We set 
n = >/2, n the refractive index, in order to study quad- 
rangle QS modes by noting the maximum openness ef- 
fects around the critical incident angle Xc = arcsin(l/n) 
[sj. In this paper, we focus on TM(transverse mag- 
netic) polarization in which wave function and its normal 
derivative are continuous across cavity boundary. Mode 
solutions are obtained by the boundary element method 
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FIG. 2: (Color online) Short-time evolution of the critical line. 
The thick and thin lines are the images of the critical line 
(green straight line) after two and four bounces, respectively. 
The good agreement with Husimi function below the critical 
line has an implication for the QS-mode formation, (a) £ — 
0.06. (b) e = 0.12. 

[Toj . Figure 1 (a) shows a QS mode found at a defor- 
mation £ = 0.12. The internal pattern looks somewhat 
complicated due to the interference between clockwise 
(red arrows) and counter-clockwise (blue arrows) waves 
propagating along each quadrangle path. The bouncing 
positions are clearly seen in the Husimi function for in- 
cident waves [20| as shown in Fig. 1 (b) where only the 
counter-clockwise part, p > 0, is shown for a clear view, 
p = sinx, X is incident angle and s is the boundary co- 
ordinate. It is evident that the bouncing positions locate 
between two periodic orbits, rectangle and diamond or- 
bits (see rectangles and diamonds), so this is a QS mode. 
The Poincare surface of section (PSOS) plotted in back- 
ground shows a chaotic ray dynamics. 

We note that the openness effects, such as G-H shift 
and FF effect, needed in the quasiscarring still exist even 
when ray dynamics is almost regular. This means that 
the quasiscarring does not demand the chaotic ray dy- 
namics and QS modes can be found in the small defor- 
mation regime with almost regular phase space. This is 
the case shown in Fig. 1 (c), a QS mode at £ = 0.06. 
The bouncing positions shown in Fig. 1 (d) clearly show 
that the corresponding quadrangle is not a periodic orbit. 
One can see almost regular ray dynamics from the PSOS 
in background. It is emphasized that the QS modes ob- 
served cannot be explained as a scarred mode modified 
by openness effect such as G-H shift and FF effect which 
give small shifts of a wavelength order and a degree in in- 
cident angle [3]. Therefore, as in the spiral-shaped cavity 
there is no periodic orbit associated with the QS mode 
pattern. 

The bouncing points of counter-clockwise waves in QS 
modes of Figs. 1 (b) and (d) locate on the left side 
of those of the diamond orbit. Why not on the right 
side? In order to understand this point, we have to take 
short-time ray dynamics into account. It is emphasized 
that even in a chaotic case rays being above the crit- 
ical line come down below the line and escape follow- 
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FIG. 3: (Golor online) (a) All existing modes (gray dots) in 
a circular case and the modes of even-even symmetry class, 
( Ai , • • • , Aio , 5i , • • • , ^8 , Ci , • • • , Cs ) . They are taken as ini- 
tial modes for the evolution study of (b). Diamonds and tri- 
angles imply the mode patterns at a larger deformation, dia- 
mond and QS mode patterns, respectively, (b) Evolution of 
Im[A;i?] of the chosen initial modes as increasing e. Existence 
of an EP at £ 0.05 is seen as a singular-like behavior of Aq 
and Be, and another EP around e ^ 0.073 (see bold lines and 
arrows), (c) Transition from crossing of (A5, B^) and (Ae, Bq) 
to avoided crossing of (A7, B7) of Re[A;i?] at the EP. (d) Mode 
patterns of a QS mode and a diamond mode, corresponding 
to green and blue dots in (b). 

ing some dynamicalpassages associated with unstable- 
manifold structure [sl, |21[ . Let us consider a ray ensem- 
ble evenly distributed above the critical line and see how 
the rays evolve in a short-time period, especially near the 
critical line where QS modes are localized. It can be il- 
lustrated by evolution of the critical line as shown in Fig. 
2. The green and yellow lines denote the images of the 
critical line after two and four bounces, respectively. It is 
clear that the passages to escape locate on the left side of 
the diamond orbit, and they are in a good agreement with 
the bouncing points of the QS quadrangles. This implies 
that the pattern of the QS modes is the consequence of a 
natural constructive interference, under openness effects, 
of the waves coming down from the upper part of the 
critical line. 

In fact, the appearance of QS modes shown in Fig. 1 
is not exceptional, but rather common in typically de- 
formed microcavities. In order to see their common exis- 
tence, we investigate evolution of the WGM modes with 
incident angles x similar to the critical angle Xc for total 
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internal reflection in the range of 52 < m < 80, m is 
the angular quantum number (m is even in our case of 
even-even symmetry class). In Fig. 3 (a) small dots are 
all WGMs existing in the range shown and each mode 
has two quantum number (m,/), / is the radial quan- 
tum number. The modes of the lowest Ke[kR]^ k the 
wavenumber, at each m belong to / = 1 mode group, and 
the next-lowest ones to / = 2 mode group and so on. The 
orange line denotes Re[kR]=m corresponding to the crit- 
ical incident angle Xc = arcsin(l/n). For convenience, 
we mark the modes of / = 4 as Ai, • • • , Aio, the / = 5 
modes as 5i, • • • , and the / = 6 modes as Ci, • • • , C3 
as indicated by diamonds and triangles in Fig. 3 (a). 

Evolutions of these modes with e are well characterized 
by the plot of lm[kR], related to decay rate, as shown 
in Fig. 3 (b). The red lines denote evolutions of Ai 
modes, and the blue and green lines do those of Bi and Q 
modes, respectively. From a low deformation about s 
0.02 almost modes begin to exhibit QS character, and we 
will call them high-Q QS modes if their incident angle is 
greater than Xc and low-Q QS modes otherwise, some of 
which are indicated by upward-pointing and downward- 
pointing triangles, respectively, in Fig. 3 (a). The lm[kR] 
values of these QS modes converge and are connected 
to two branches after e ~ 0.08, the modes in the upper 
branch exhibit a QS mode pattern and others in the lower 
branch do a diamond mode pattern. Figure 3 (d) shows 
these patterns corresponding to C2 mode and Bq mode 
at £ = 0.12 (green dot and blue dot in Fig. 3 (b)). 

Another interesting feature is the existence of two ex- 
ceptional points (EPs) near e 0.05 and e 0.073. The 
EP is a degenerate point in a non-Hermitian system like 
an open quantum system, at which two eigenvalues and 
eigenfunctions coalesce [10| . Mathematically, the EP is a 
square-root branch point in a parameter space [l^. To 
find an EP, we need two parameters, i.e., its codimen- 
sion is two. In our case we have quasi- two dimensional 
parameter space, one is the deformation e and the other 
is a discrete internal parameter m H, [22 1. Due to the 
discreteness of m, we cannot find the exact EP position, 
but can recognize the vicinity of EP from the singular- 
like behavior of the evolutions. The evolution lines of 
Aq and Bq modes show abrupt changes at £ 0.05, and 
those of A2 and C2 modes do at £ 0.073 (see bold lines 
in Fig. 3 (b)). These EPs can be also confirmed by the 
transition behavior from crossing to avoided crossing in 
Ke[kR] plot, as shown in Fig. 3 (c), as decreasing m. 
We also see that kR values of Aq and Bq are very close 
at 5 0.05. The same behavior appears near e 0.073 
in the A2 and C2 evolutions. In fact, by introducing 
another continuous parameter like the refractive index 
n, we can find nearby EPs, in our case, EPs locate at 
{riEP^eEp) ^ (1.4164,0.050) and (1.4160,0.0730). 

Due to the nontrivial topology of energy surfaces near 
an EP, neighboring two modes showing almost identical 
mode properties can be completely different modes as 
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FIG. 4: (Color online) Change of peak positions of Husimi 
functions for two interacting QS modes, the high-Q QS mode 
(red upward-pointing triangle) and the low-Q QS mode (blue 
downward-pointing triangle) before and after e = Sep- The 
ray return image of the high-Q QS mode point marked by 
green circle is shown, (a) Sep — 0-05 case of Aq and Bq 
modes, (b) Sep — 0.073 case of A2 and C2 modes. The 
critical lines are denoted by orange lines. 



a parameter passes by an EP [lO|, [ij]. This branching 
property is observed in Fig. 3 (b). For example, the EP 
existing at eep — 0.05 divides similar modes, (^6,^7) 
or (^6, ^7)5 dit £ < Sep into different modes, a QS mode 
(^6,^7) and a diamond mode (^47,^6), at £ > Sep- As 
a result, while (A3, ■ ■ ■ ^ Aq^ B^^ B^) modes evolve into the 
QS branch, (^7, • • • , Aio, Bi^ ■ ■ ■ ^ Bq) modes do into the 
diamond branch, as marked in Fig. 3 (a) by triangles and 
diamonds. Similar behavior is expected to happen by the 
other EP at £ep — 0.73, while in Fig. 3 (b) the modes 
passing by the EP look approaching the same QS mode 
branch. Actually, looking into details, there are two kinds 
of QS modes at e > eEp{— 0.73). They have a small 
difference in bouncing positions and Q factor, so we say 
QS-I modes and QS-II modes for convenience. The QS-I 
mode at £ = 0.12 has the pattern shown in Fig. 1(a), and 
the other QS-II mode has slightly different pattern shown 
in the upper panel of Fig. 3 (d). Therefore, {Ai^ A2^Cz) 
and (A3, • • • , Ae, Ci, C2) become the QS-I and the QS-II 
modes, respectively, as marked by red and blue colors in 
Fig. 3 (a). 

The interaction between a high-Q and a low-Q QS 
modes can be modeled by a two-by-two interaction 
Hamiltonian matrix with diagonal elements E^^h.l] = 
^{H,L} ~ ^7{i^,L}/2, (jO is the frequency and 7 the decay 
rate, and the off-diagonal element Chl representing cou- 
pling strength between them [l^]. The condition for an 
EP is then given by Chl = — 7l|/4 when uoh = ^l- 
In our case, all elements change depending on e. For a 
small deformation, \'^h — 7l| is much greater than AChl^ 
similar to an integrable case. As increasing deformation, 
not only \'^h — 7l| decreases (see Fig. 3 (b)) but also 
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FIG. 5: (Color online) Mode evolutions in the stadium de- 
formation. All modes can be classified into three modes, a 
diamond scar mode and two QS modes whose patterns are 
shown at the right panels. 

the Chl increases due to introduction of nonintegrabil- 
ity of ray dynamics. As a result, at some deformation 
the condition Chl = — 7l|/4 would be hold, and the 
point would be an EP if uoh = cjl at the deformation. 
Therefore, the EP can exist roughly when the coupling 
Chl increases drastically. 

The change of Chl can be qualitatively estimated from 
the ray return image of the high-Q QS pattern. For 
example, consider Aq mode at e = 0.02 whose a peak 
position of Husimi function is the red upward-pointing 
triangle in Fig. 4 (a). The rays at the position (see the 
green circle) return, after four bounces, to the left side 
of the initial position as indicated by the green image, 
implying a small coupling Chl to the low-Q mode {Bq 
mode) existing below (see the blue downward-pointing 
triangle). At a higher deformation e = 0.04, the return 
image locates below the initial point due to the widen 
nonlinear resonance. This indicates more larger coupling 
with the low-Q QS mode since wave functions of Aq and 
Bq modes are expected to show a larger overlap. This in- 
crease of Chl is consistent with the existence of an EP at 
e c± 0.05. After the EP, e.g., at £ = 0.06, the peak points, 
initially showing a vertical difference (a large Q differ- 
ence), distribute with a horizontal difference (a small Q 
difference), the high-Q QS mode becoming a QS-II mode 
and the low-Q QS mode approaching a diamond mode. 
The peak variations and the return image of A2 and C2 
modes exhibit the same behavior, as shown in Fig. 4 
(b), except for the different final modes, QS-I and QS-II 
modes. 

The branching role of the EPs for QS modes are quite 
generic. As shown in Fig. 5, a similar evolution behav- 
ior of lm[kR] can be found in the stadium deformation 
showing a fully chaotic ray dynamics. The deformation 
parameter 77 is defined by 77 = a/2R^ a is the length 
of the parallel boundary and R is the radius of circular 
boundary. In this stadium case, A^ and ^5 mode evolve 



through the vicinity of an EP located around r] = 0.025 
associated with two branches, a diamond scar branch and 
QS-II branch, as shown in Fig. 5. A3 and C3 modes also 
pass by another EP near 77 = 0.06 related to two QS 
branches, QS-II and QS-I branches. The modes patterns 
of the three branches are shown in Fig. 5. This branch- 
ing behavior by two EPs is similar to the quadrupole 
case, i.e., (^1,^2,03,04), (A3, A4, ^5, • • • , ^8, Oi, C2), 
and (As, • • • , Aio, 5i, • • • , B^) modes arrive at the QS-I, 
the QS-II, and the diamond branches, respectively. Note 
that the modes in the diamond branch have higher Q 
values, unlike the quadrupole case, due to the different 
boundary geometry, and show small wiggles implying in- 
teractions with other high-Q modes. This similar branch- 
ing behavior of the stadium case means that the inter- 
action between QS modes is not directly relevant to the 
PSOS structure. Instead, the combination of openness 
and short-time ray dynamics is more responsible for the 
interaction. Other branching behavior by an EP is ex- 
pected for other modes since it is based on the singular 
feature of the EP. 

In conclusion, we have shown that appearance of QS 
modes is common in a typically deformed microcavity 
such as a quadrupole and a stadium microcavities. At a 
small deformation, they can be divided into high-Q and 
low-Q QS modes, and both QS modes evolve into three 
branches, a diamond branch and two QS (QS-I, QS-II) 
branches. This branching behavior has been understood 
by the existence of two EPs, i.e., the inherent property 
of the square-root branch points. The mode pattern of a 
QS mode can be explained by a combination of openness 
and short- time ray dynamics near the critical line, rather 
than its PSOS structure, due to its leaky nature. 

This work was supported by NRF grant (2010- 
0008669). 
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